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ABSTRACT 

We investigate the prospects of detecting radio afterglows from long Gamma-Ray 
Bursts (GRBs) from Population III (Pop III) progenitors using the SKA precur¬ 
sor instruments WMA (Murchison Widefield Array) and ASKAP (Australian SKA 
Pathfinder). We derive a realistic model of GRB afterglows that encompasses the 
widest range of plausible physical parameters and observation angles. We define the 
best case scenario of Pop III GRB energy and redsliift distributions. Using probability 
distribution functions fitted to the observed microphysical parameters of long GRBs, 
we simulate a large number of Pop III GRB afterglows to find the global probability of 
detection. We find that ASKAP may be able to detect 35% of Pop III GRB afterglows 
in the optimistic case, and 27% in the pessimistic case. A negligible number will be 
detectable by MWA in either case. Detections per image for ASKAP, found by incor¬ 
porating intrinsic rates with detectable timescales, are as high as ~ 6000 and as low 
as ~ 11, which shows the optimistic case is unrealistic. We track how the afterglow 
flux density changes over various time intervals and find that, because of their very 
slow variability, the cadence for blind searches of these afterglows should be as long as 
possible. We also find Pop III GRBs at high redshift have radio afterglow lightcurves 
that are indistinguishable from those of regular long GRBs in the more local universe. 
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1 INTRODUCTION 

Gamma-Ray Bursts (GRBs) are the brightest explosions in 
the universe. These transient events outshine entire galax¬ 
ies for a short time. The degree of variability and short life 
of GRBs points towards a compact source. Typical isotropic 
equivalent energies lie in the region of 10 53 erg. Such high en¬ 
ergy from a compact object can be explained if th e emission 


is narrowly collimated and highly relativistic ( 

Rhoad^ 1997|; 

Sari et al. 1999; Rhoads 1999; Frail et al. 2001 

: Bloom et al. 

20031: Chandra & Frai| 2012j). The most plausible hvpothe- 


sis for the production of a GRB is a system consisting of 
a stellar mass black hole surrounded by an accretion torus 
within a magnetic field. Angular momentum in the accre- 
tion torus is extracted via t he Blandford-Znajek mechanism 
iBlandford fc Znaiekl 1 19771 ) to produce a polar outflow of 
highly relativistic plasma. Velocity variations in the outflow 
lead to internal shocks where the electrons are accelerated 
to produce gamma-ray photons. 

In addition, the outflow produces external shocks when 
it interacts with an ambient medium, leading to synchrotron 
radiation (afterglow of the GRB). The peak frequency of 


this afterglow decreases as the outflow loses energy; typi¬ 
cally what begins as a predominantly X-ray afterglow, over 
time becomes mostly optical, infrared, and eventually a ra- 
dio a fterglow (iMeszAros fe Reesll 19971 : ISari et abll 1 9981 : iPiranl 
120041 : IZhand 120071 : Gao et ahl 20131 ) . A similar mechanism 
explains the relativistic outflow and radio lobes seen in ac¬ 
tive galactic nuclei (AGN), with the difference being mainly 
one of scale: AGN involve supermassive black holes, are less 
bright, less variable, and much longer lived dMirabell 2004]; 

erg et all: 


lGhisellinll2005l : iNemmen et al.ll2012r iBromberg et al .1120 1 il l . 

GRBs are classified as either short/hard or long/soft 
(as well as more exotic types such as under-luminous and 
ultra-long). Within the framework of GRBs being produced 
by black hole formation/accretion, the distinction between 
the two categories can be explained by fundamentally dif¬ 
ferent methods of formation or accretion onto a stellar mass 
black hole. Short duration GRBs with hard gamma-ray 
spectra are thought to originate from a merger of com¬ 
pact objec ts, two neutron stars or a black hole and a neu¬ 
tron s tar (iGoodman 19861 : Woosle^ 19931: iRuffert fe Janka 


119991 : iBloom et al.l 20021 : lAlov et al. 2005l ~ Rezzolla et al 
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1201 ll'l . Long duration GRBs with soft gamma-ray spectra 


tivcly massive star (IWooslev 1993}: 

MacFadven & Wooslev 

19991: Frail et al. 2001; Bloom et al. 

2002|). 


Short GRBs are unsuited as probes of Pop III stars, as 
the production of the compact stellar remnant progenitors 
introduces a large degree of separation between the Pop III 
star and the short GRB. Where one progenitor is a black 
hole, the only attribute connected to the original star is its 
mass. Add to this the time taken to produce the compact 
remnants, but mostly in inspiral time for the two remnants 
to merge, short GRBs are too far removed fr om the origi- 
nal stars to be suitable probes of Pop III stars (IBloom et al.l 
l2002fh Furthermore, short GRBs tend to have very faint af¬ 
terglows. 

On the other hand, the rapid evolution and short life 
of a massive star which culminates in a long GRB links 
those GRBs closely with the initial star formation. Therefore 
long GRBs from Pop III stars can be used to constrain the 
formation history of the first generation of stars. 

The potential of utilising GRBs as tools to study Pop 
III stars has been strengthened by observ ations of ve ry hig h- 
z afterglows (eg GRB090429B z ~ 9.4 (ICucchiaral 1201 IT ) I . 
The formation history of Pop III stars has been theoretically 
studied by multiple authors. Such work includes the ionisa¬ 
tion rate of the IGM and universal chemical evolution his¬ 
tory that explain star formation rates in the local universe. 
Theoretical Pop III star formation rates (SFRs) peak be- 


2001} 

; 

Maio et al. 

2012 

; 

Wise et all 


120101 ; Ide Souza et al.l 1201 ll : I Johnson et al.l 


rate at which stars form out of the primordial environment, 
with the e arliest limits of Pop III sta r formation being placed 
at z ~ 60 (INaoz fe Bromberg! l2007l h In the most optimistic 
case the peak of formation is at z ~ 6, with appreciable for¬ 
mation still occurring at z ~ 3. With GRBs detectable at 
z ~ 9.4, it should be possible to use GRBs to probe what 
is possibly the peak of Pop III star formation, or at least 
to constrain the redshift limits at which Pop III formation 
ceased. 

As a foundation for our simulation, we employ previ¬ 
ous studies to derive a model that is capable of accepting 
the broad distribution of physical parameters that determine 
the afterglow lightcurve. We refine the method of account¬ 
ing for spherical shape and beaming of the afterglow emis¬ 
sion region to improve the efficiency of the calculations, and 
expanded it to allow for arbitrary orientation angles. Our 
refinements result in a new model that can take practically 
any set of plausible physical parameters of a GRB, at any 
orientation angle and observing frequency, and return the 
lightcurve over 10 12 seconds within a minute on a desktop 
computer. 

We have adopted this approach (see Macpherson et al. 
2013) in favour of the construction of a luminosity func¬ 
tion for Pop III GRB afterglows for three reasons. Firstly, 
given that Pop III GRBs are still highly speculative objects, 
constructing a realistic luminosity function is ambiguous. 
Secondly, the luminosity function considers the brightness 
at a specific time, usually one day after the burst. Most ra- 
dio afterglows are still increasing in brightness at this time 
(IChandra fe Frail 12012 1. thus using a luminosity function, 
one would underestimate the proportion of detectable af¬ 


terglows. Thirdly, the suitability of luminosity functions to 
transients like GRB afterglows is dubious; they are more 
suited to objects whose luminosity can be considered con¬ 
stant or stable for most intents and purposes. 

The paper is organised as follows: the afterglow simula¬ 
tion method is described in Section 2, and Section 3 de¬ 
scribes the distributions of GRB parameters used in the 
model. The results of the simulations of Pop III GRB af¬ 
terglows are presented and discussed in Section 4. Section 
5 summarises our main results, the detectability of Pop III 
afterglows by SKA precursors. 


2 AFTERGLOW SIMULATION METHOD 

As a foundation to simulate the light curves of the GRB 
afterglows we use a method based primarily on Mesler et al. 
(2014). We begin by calculating the initial parameters of the 
ejecta: the bulk Lorentz factor (Po), the ejecta mass (M e j), 
and initial radius ro. 

The formulation we use to calcul ate th e initial bulk 
Lorentz factor To is taken from iToma et al.l (l201lli . with 
an initial time to set to 1 second, 


r 0 


£ iE 


\4nnm p c 5 tQ 


1/8 


(1) 


where Ei BO is the isotropic equivalent energy of the GRB, m p 
is the proton mass, and n is the external medium density. 

M e j, and the deceleration radius r 0 at which T be- 
gins to evolve, are found using the method described in 
IPanaitescu fe Kumail d2000h : 


( (3 -s)E \ 1/3 

\4:TvAmpC 2 r 2 J 


( 2 ) 

(3) 


In the above equations, E is the total energy of the 
fireball, which is found by scaling Ei so by the fraction of a 
hemisphere projected by the initial jet opening angle 9j , that 
is E = E iso (l — cos 9j). The values A and s in equation [3] are 
properties of the external medium density, i.e. n(r ) = Ar~ s . 
In the case of an isotropic medium, s = 0 and A = n. 

The ma ss mp swept up by the fireball by radius rp 
is fou nd by llHuang et al.lll999l ; IPanaitescu fe Kumaill2000l : 
|Pe ; erll2012h : 


m{r 0 ) = (4) 

The evolution of F is evaluated numerically over inter¬ 
vals of r, for which the rate of change of V is: 
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dr 

dm 
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t\ 

B 


syn 

/2 


7 (r 2 -1) - (7 - i)r^ 2 


M e j + em + (1 — e)m(27r — (7 — 1)(1 + T~ 2 )) 

4r +1 


3r 




e t , ~ 1 + 

°syn 1 ^exp 

67r m e c 


ot-B' 2 7 m ’ 

327reflwn p c 2 r(r — 1), 


(5) 

( 6 ) 

(7) 

(8) 
(9) 


where 7 is the adiabatic index of the jet material, t is the 
radiative efficiency, e e is the fraction of blast wave energy 
in the ele ctrons, f' v „ is the comoving s ynchrotron cooling 
timescale dDai et al. Hl999l: fRh oa dsl ll999lh texp is the comov¬ 
ing age of the jet, ctt is the Th ompson cross se ction, B' is 
the comoving magnetic field IlMesIer et al.|[2014ll . and eg is 
the magnetic field energy fraction of the blast wave. 

The minimum Lorentz factor of injected electrons, 7 m , 
is calculated in different ways depending on the value of the 
electron energy index, p. Where p>2, 


^fm — 


l + e e (r-l) 


m p (p - 2 ) 
m e (p- 1) ’ 


( 10 ) 


djohannesson et al.l 120061 : iMesler et al.l l2014h . If the 
electron energy distribution is unusually flat (p< 2 ), 


7m = 1 + 


(e e (r-i) 7 r 2 


m p (2 - p) 
m e (p - 1 ) 


7a - 



i/(p-i) 


( 11 ) 

( 12 ) 


Equation |TT| dBhattacharval 1200 ll : IDai &: Chend 1200 il l 
depends upon the maximum injected electron L orentz facto r 
7 u, here approximated by Equation 1121 dToma et al.ll 20 lHl . 
where e c is the electron charge. 

The differential mass swept up dm is approximated by: 


—— = 4nr 2 n(r)m p (13) 

dr 

dHuang et al.lfl999l : IPe’eill 2012 h . With these eq uations 
we find T and m in terms of r, which IMesler et all d2014l l 
relate to local rest frame time t as: 


dr c/ 3r(r + vr 2 - 1)' 

The evolving opening half angle 9j is found by calculat¬ 
ing Csj the_sgeed of sound within the blast wave (Equation 
1151 IMesler et aid d2014h l. Over a time interval dt, the blast 
wave spreads sideways by an amount c s dt., in addition to ex¬ 
pansion due to radial motion. Defining a as the lateral size 
of the blast wave, we calculate 9j after a small change in 
radius dr by: 

© 0000 RAS, MNRAS 000, 000-000 


c s = cT 


7 ( 7 - i)(r- 1 ) 

1 + 7(T-1) 


a = rtan( 0 j) + (c s dt), 
9j = arctan(a/r). 


(15) 

(16) 
(17) 


The synchrotron emission spectrum of the blast wave 
depends on three characteristic Lorentz factors: the mini¬ 
mum injection Lorentz factor y m , the cooling Lorentz factor 
7 C , and the absorption Lorentz factor 7 ^. 

The Lorentz factor of cooling electrons is 


6nm e c 


B' 2 a T t'(Y + l) 


(18) 


|Sari_et_aL _ 1 19981: _ iBhattacharval 


2001 


Johannesson et all 120061 : iToma et ~ I l201ll : IMesler et aid 


2014lh where Y is the Compton parameter, whose value 


depends on the relative values of the characteristic Lorentz 
factors as shown below. 


{ 7m7cT-e 7m > 7c 

T- e 7^T 1 7c 3_P) 7c > 7m 

r e (C'2 _P 7l7m P_1) ) 1/(p+5) 7a > rnax( 7 m , 7 c), 

where r e is the optical depth: 

(jTm(r) 


(19) 


47 t m p r' 2 ’ 


and 


Ci = 


5 e c T e 
otB' 


( 20 ) 


( 21 ) 


Because the equations of the Compton parameter Y de¬ 
pend on the relative values of 7 C , 7 m , and 7 a , one must find 
the solution to 7 C which is self-consistent. This is compli¬ 
cated further by the fact that also depends on the rel¬ 
ative values of the other characteristic Lorentz factors, as 
described below: 


7a = 


,-<0.3 -1/2 
02 7c 

C 2 03 7m 1/2 

(C 27 c) 1/6 

(C27)T 1 ) 1/(p+4) 

(Cl27c7)r 1 ) 1/(p+5, 


7a < 7c < 7rn 
7 a ^ 7 m ^ 7 c 
7c < 7a < 7m 
7m < 7a < 7c 
7a > max(7m,7c), 


( 22 ) 


We approach this problem by finding the values of y c 
in all possible scenarios, then select whichever one is self- 
consistent. For instance, if the equation for Y in the fast 
cooling regime (i.e. > 7c) results in 7 C greater than 7 m , 

obviously this solution is invalid and the jet is not in the fast 
cooling regime. Then, if the Y equation for the slow cooling 
regime also yields > 7m, the jet is definitely in the slow 
cooling regime and the value of 7 c is valid. 

Similarly in finding the 7 a , we find all possible values, 
under every arrangement of the characteristic Lorentz fac¬ 
tors, and find that which is self-consistent. 

These characteristic Lorentz factors transform to char¬ 
acteristic frequencies via: 
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j/( 7 ) 


3Te c B' 

(1 — /5)47rm e c’ 


(23) 


where /3 is the velocity (in terms of c) corresponding to the 
Lorentz factor 7 . 

Then the peak flux density is calculated via 


(1 + z)y/3<j> p el 
d-Kd 2 L m e m p c 2 




(24) 


where <j> v is t he p -dependent factor defined in 
IWiiers fc Galamal (119991 . dz, is the luminosity distance 
to the GRB, and /3 m is the velocity of the lowest-energy 
injected electrons. Prior to the jet break at time tj, F™ ax 
remains approximately constant. After the jet break, the 
peak flux den sity evolves wit h time as J 7 ’™ a ' x — /t) 

j|Rhoadslll999l : lBhattacharvall200l| j . where is the peak 
flux density calculated at the point of the jet break. 

The jet break marks the stage at which the centre of the 
jet head becomes causally linked with its outer edge; when 
the bulk Lorentz factor of the jet is less than the inverse of 
its angular separation (r < 1 /Oj). In the model we use, 9j 
evolves with all other quantities. This evolution is slow at 
early times, accelerating at the point of jet break. 

The peak flux density is emitted at the median of the 
three characteristic frequencies. The slope of the spectrum 
between characteristic frequencies depends upon those fre¬ 
quencies’ relative values, as detailed below (all frequencies 
in the frame of the local ISM). 

I 11 the fast cooling ( v c < v m ) regime: 


F v = F„ max 


(v/v a ) 2 (v a /v c ) 1/3 

(Cvc) 1/3 

Cc/v) 1/2 

{vlVm)~ V/2 {v c /Vm) 1/2 

C/VcfCc/Va .?' 2 

C/^f 2 

K/M 1/2 

HVn,)- P,2 CalVm) 1/2 
C v/v c ) 2 (v c /v a ) 3/2 
iy/vaf' 2 

(y/va)~ p/2 


V < V a < V c < V m 
V a < V < V c < V m 
V a < V c < V < V m 
V a < V c < V rn < V 

V <C V c < V a Vm 
v c ^ <C V a Vm 
V c < V a < V < V m 
V c V a <1 Vjn <C V 

V < V c < V m < V a 

V c < V <V a 
Vc < V m < V a < V 

(25) 


In the slow cooling ( v m < v c ) regime: 


(v/v a ) 2 {Va/Vm) 1/3 

V < V a < V m < Vc 

{v/Vm) 1/3 

V a < V < Vm < V c 

(: v/v m )- ( p-D /2 

Va < V m < V < V c 

(V/Vc)- P/2 (VC/Vra)- (P - 1)/2 

Va < V m < V c < V 

{V/V m ) 2 {v m /V a f /2 

V < Vm < V a < V c 

(CVa) 6 ' 2 

Vm < V < V a < V c 

(v/Va) (P ~ 1)/2 

Vm < V a < V < V c 

(v/Vc) P/2 (Vc/Va) (P ~ 1)/2 

Vm < V a < V c < V 

{v/Vm) 2 (Vm/Va) 5/2 

V < V m < V c < V a 

(CVa.) 312 

Vm < V < V a 

(v/v a )~ P/2 

Vm < V c < V a < V 


(26) 

Hence the flux density is obtained at any given fre¬ 
quency v for a local observer. However, we have not yet 


accounted for the spherical nature of the emitting region 
or its relativistic velocity. The spherical shape introduces a 
slight delay to radiation emitted from angles further from 
the observer’s line of sight. Were the flow not relativistic, 
this delay could be calculated from simple geometric argu¬ 
ments; effects of special relativity complicate this angle de¬ 
pendent delay. Further, the Doppler beaming factor for high 
T amplifies emissions along the line of sight while reducing 

emissions at angles greate r than 1/r. _ 

From the appendix of lMesler et alj J2014h the equal ar¬ 
rival time surface is 


f (1 — /3cos$) 

t = constant = / --- dr, 

J c/3 


(27) 


where it is suggested that this equation be used to find the 
emission radius, r, of radiation from a section of the jet at an 
angle 9, as measured from the explosion centre, to the line 
of sight, by integrating over an increasing r with constant 
9 until reaching the observation time t. This must be done 
for 0 < 9 < 9j to find the equal time arrival surface r as a 
function of 9. 

Repeated integration of equation ([27|) to find r for 
incremental values of 9 is a time consuming process. It 
is more efficient to do the opposite: to find 9 for incre¬ 
mental values of r. The equation can be rearranged into 
t = —rcos9/c + f dr/c/3, and from there it can be shown 
that 



(28) 


Thus radiation emitted at some smaller radius n = r — 
Ar, reaching the observer simultaneously with that emitted 
at radius r, will have been at an angle 9 to the line of sight. 
The total observed flux density is given by 


F v 


fT 


d cos 9 d(j>, 


(29) 


where L' v ,[r(9)\ is the total radiated luminosity at the co¬ 
moving frequency v' at r(9), L' v , = Fu/T, T> is the Doppler 
factor T> = 1/T(1 — /3 cos 9), and Q.j is the jet solid angle. 

When the observer is oriented precisely on the jet axis, 
f d(j) = 2n, as the emitting region appears a full circle. In 
more general terms, if the observer is offset from the jet axis 
by some angle 9i, f A(j> = 2n only for angles 9 ^ 9j — 9i, 
and evaluates to 0 when 9 < 9i — 9j and 9 > 9i + 9j. When 
9i — 9j < 9 < 9i + 9j, f dcj> = 2arccos ((9 2 + 9 2 — 9 2 )/299i). 

Thus we iteratively reduce r by an appropriately small 
fraction Ar, for each value calculating the emission angle 9. 
As long as 9 < 9i + 9j , we also calculate the total luminosity 
= ^( r )/r(^), Doppler factor V(r,9) = 1/F(r)(l — 
/3(r)cos0), and f d</>. Once done, we integrate L' v ,[r[9)\D 3 
from 0 to 9i 9j , and divide by fito find the total observed 
flux density. 

Figure |T] shows how the model’s afterglow lightcurves 
compare with observations in both a constant density 
medium and a wind medium. The common parameters of 
the simulated GRBs are: Ei SO = 10 53 erg, z = 1, 9j = 0.1 
rad, p = 2.3, e e = 0.1, eb = 0.01. In the case of the isotropic 
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Figure 1 . Comparison of observed (dotted lines) and simulated 
fiducial GRB 8.5GHz radio afterglow lightcurves in the case of 
an isotropic external medium (solid line) and a wind medium 
(dashed line). 


medium, the density n is lent -3 . In the case of the wind 
medium, the density is inversely proportional to r 2 , and the 
values of the deceleration radius ro and density at that ra¬ 
dius n(r o) are set to the fiducial values in Panaitescu & 
Kumar (2000, equations (8) and (9)) for a GRB from a Wolf- 
Rayet star. These are plotted over a sample of observed GRB 
radio afterglow lightcurves at frequencies at or near 8.5GHz. 
The data for these radio afterglows was provided by Poonam 
Chandra (2014, private communication). 

The simulated lightcurves in Figure |Tj are comparable 
to observations, and this shows that the simulation is real¬ 
istic. An Eiso of 10 53 ergs and redshift of 1 applied to both 
cases (constant density and wind medium) are roughly me¬ 
dian values as measured through observations, and produces 
roughly median lightcurves. We do not consider whether this 
Ei so and z are true medians of the intrinsic distributions, as 
both the lightcurves and the apparent parameters are drawn 
from observations, hence are subject to the same selection 
effects. 

The differences between the wind density and constant 
density medium cases are small. The lightcurve in the wind 
density case has a steeper rise to peak earlier than in the 
constant density case. At later times, and for the majority of 
the time when the afterglows are observable, the lightcurves 
tend to follow each other with less than half an order of 
magnitude between their flux densities. 

We will only consider constant density media in the re¬ 
mainder of this paper. To properly construct a wind density 
profile requires a robust estimate of the mass loss rate of 
the progenitor star. In the case of Pop III stars, this is far 
from trivial. It has been argued that the rate of radiation- 
driven mass loss from a Pop III star should be minimal due 
to the lack of absorption lines in its mostly hydrogen enve¬ 
lope and environment. On the other hand, it is argued that 
in order for a massive Pop III star to collapse while retaining 
enough angular momentum to form a GRB, it must undergo 
chemically homogeneous evolution (jYoon et al.ll2012l ), which 
would dredge heavier elements produced in the core up to 
the outer layers, restoring radiation driven mass loss. Fur¬ 


thermore, chemically homogeneous evolution requires the 
star to have a high rotation velocity, which implies that the 
kinematic mass loss should be high. This last point may or 
may not be influential to GRB afterglows, as the kinemat¬ 
ics that drive mass loss are greatest at the star’s equator, 
and zero at the poles where the relativistic GRB jets are 
expected to launch. 

Finding a satisfactory solution to the problem of the 
properties of a wind medium around a Pop III star would 
be a separate and substantial investigation. Given the rela¬ 
tively small difference in the resultant lightcurves as shown 
in Figure [I] we keep to the simpler case of assuming constant 
density media. 


3 DISTRIBUTIONS OF GRB PHYSICAL 

PARAMETERS 

The relevant physical parameters that determine GRB after¬ 
glows are the total isotropic equivalent energy f?i so , redshift 
z, the jet half-opening angle 9j, density of the surround¬ 
ing medium n, the electron energy distribution index p, the 
electron energy fraction e e , and the magnetic field energy 
fraction es. 

We have no reason to expect that the distributions of 
the parameters 9j, n, p, e e , or es for Pop III progenitor 
GRBs will differ from those of normal GRBs. Many GRB 
simulations use a set of ‘fiducial’ parameters as a way of esti¬ 
mating typical GRB lightcurves. In this work we investigate 
GRBs where the parameters are allowed to vary as observa¬ 
tions show, thus creating a realistic set of GRB lightcurves. 
The fiducial parameters provide a starting point and reality 
check for the param eter d istrib uti ons. A typical set of fidu¬ 
cial parameters from 1 Toma et al] (l201ll ) is 9j = 0.1, n = 1.0 
cnF 3 , p = 2.3, e e = 0.1, and es = 0.01. 

We do expect that the distributions of Ei so and z are 
closely linked to the distribution of progenitor mass and the 
SFR, respectively. As these distributions are uncertain, we 
investigate two extreme scenarios for the Ei ao and z dis¬ 
tributions: under the most optimistic and most pessimistic 
combinations of Pop III initial mass function (IMF) and 
SFR. 


3.1 GRB isotropic equivalent energy Ei S 


The energy reservoir of a long GRB is associated with the 
mass of its progenitor star. The amount of energy released 
is proportional to the mass of the accretion torus around 
the nascent black hole, which in turn is pr o portional to the 
total mass of the star d Toma e t al.ll20lH ). [Nakauch i et al.l 
(120121) and lSuwa fe Iokal (1201 ll) simulated the production of 
GRBs from collapsing Pop III stars of various initial masses, 
and found the resultant Ei S0 s. The results were summarised 
in Macpherson et al. (2013; Figure 1), from which a linear 
relation between progenitor mass and resultant GRB Ei so 
can be seen, described as 


Eiso(x 10 55 erg) = O.O124(M/M 0 ) + 0.6507. (30) 

This allows us to relate the energy distribution of GRBs 
to the IMF of the progenitor stars. 

However, we do not expect every Pop III star to produce 
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Figure 2. Normalised F) so distributions of Population III GRBs 
under two scenarios of t he Pop III IMF; a high mean gaussian 
ijScannapieco et al.ll2003l l (solid line) and the Salpeter power law 
(dashed line). 


Figure 3. Pop III GRB redshift probability density functions 
under the most optimistic (solid line) and the most pessimistic 
(dashed line) Pop III SFRs. 


a GRB. According to the work of lYoon et al.l d2012l) . there is 
a certain window in mass and rotation speed in which GRB 
formation is possible. The progenitor must have an initial 
mass between ~ 12 - 84 Mq , with a r otation velocity clo se to 
the break-up limit. From figure 12 of lYoon et al.l d2012l ) and 
assuming a mass-independent distribution of Pop III stellar 
rotation speed, we generate a function of the probability of 
a Pop III star being able to produce a GRB as a function of 
mass. We find that the probability of a Pop III star forming 
a GRB is highest for a mass of ~ 60Mq. 

The IMF of Pop III stars is unknown. Previous stud¬ 
ies provide several possible alternatives to the standard 
Salpeter power la w. These include lognormal distributions 
( Tumlinsonl 20061) an d high mean Gaussian distributions 
( de Souza et al.1 1201 ll ). Of those we found it was the Gaus- 
sian distribution with / x = 55M(7) and a = 15Mq o f 
IScannapieco et aD ((2003) as used in Ide Souza et al.l d201ll ) 
that was the mo st top heavy and, combined with the 
lYoon et al.l d2012l ) mass and rotation criteria, provided the 
greatest overall probability of Pop III GRB formation. 

Assuming a mass-independen t stellar rota t ion d istribu- 
tion, we derive from figure 12 of lYoon et al.l l|2012l ) mass- 
dependent rotation limits for GRB formation. Combining 
this with the IMF and using the progenitor mass--Ei SO rela¬ 
tion, we generate the Pop III GRB Ei so probability function. 
A graphical representation of the normalised Pop III GRB 
Ei a o distributions is shown in Figure [2] 


where ’ll* (z) is the Pop III SFR, and cosmological parame¬ 
ters f lm and 12 a are 0.267 and 0.734 respectively. 

The highest Pop II I SFR fou nd in the literature is the 
most optimis tic case of de Souza et al.l (1201 if) . while the low¬ 
est is that of lKulkarni et alT i 20131) .The normalised Pop III 
GRB redshift distributions based on these SFRs are shown 
in Figure [3] 


3.3 Initial jet opening angle 9j 

Analysis of the observed opening angles of GRBs has sug- 
gest ed that #,■ evolves, wher e angles are narrower at high red¬ 
shift (%netoku3S3iI] (|2005l ). However, there is no clear phys¬ 
ical reason for this evolution, and it is alternativel y possible 
that the apparent evolution is a selection effect (iLu et al.l 
I2OI2I) . GRBs with narrower opening angles have a higher 
apparent brightness, and are therefore visible from greater 
distances. At a high redshift, an observer is less likely to 
see a GRB with a larger opening angle, so the population 
of GRBs at tha t redshift woul d appear dominated by those 
with smaller 9j. ILu et all (120121 ) found using a simulated set 
of GRBs that the apparent redshift dependent 9j distribu¬ 
tion could be reproduced if the intrinsic 9j distribution were 
lognormal, with median e M = 0.0537 and standard deviation 
a = 0 . 6 . 

This median 9j is slightly lower than the fiducial value, 
but not by an egregious amount. We therefore adopt this 
distribution in our simulations. 


3.2 Redshift distribution 

Progenitors of long GRBs are massive stars, whose lives are 
relatively short. The occurrence of long GRBs is therefore 
closely associated with the SFR. We can use theoretical es¬ 
timates of the Pop III SFR to find the redshift distribu- 
tion of Pop III GRBs. Based on the approach described in 
iBromm fc Loebl (l2002h . the non-normalised GRB redshift 
probability function is: 


3.4 Circumburst density n 

From Panaitescu fc Kumar! ll2002l); Chevalier fc Lil ll2004h : 
iPiro et al.l ( 20051 ): ICenko et al.l d 201 ll ) we take circumburst 
density measurements of twelve GRBs, shown in Table [I] 
Fitting to a lognormal distributiorQ yields distribution 
parameters e^ = 0.88 and a = 3.5. The Kolmogorov-Smirnov 
(K-S) test of the data to this distributiorQ shows a 90% 
probability that the data came from this distribution. The 


-Pgrb^) 


(1 + z) 3 \/f2m(l + Z) 3 + I2a 


d z, 


(31) 


1 Python 2.7.3, scipy 0.14.0 
scipy.stats.lognorm.fit(data, Hoc = 0) 

2 scipy.stats.kstest(data, ‘lognorm’, args = (3.5, 0, 0.88)) 
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Table 1. Observationally derived microphysical parameters of 
GRBs 


GRB Density (cm 3 ) 

P 

e e 

es 

970228 


2.44 8 



970508 

0.75 2 

2.29 2 . 3 .8 

0 . 19 1 ’ 2 ’ 3 

0.1685 1,2 ’ 3,6 

971214 


2.2® 



980329 


2.69 3,8 

0 . 12 3 

0.17 3 

980519 

0.14 2 

2.87 2 ’ 8 

0 . 11 2 

0.000035 2 

980703 


2.64 3,8 

0.27 3 

0.0018 3 

990123 

0.0019 4 

2.135 2,8 

0.13 2 

0.00074 2 

990510 

0.29 4 

1.945 2,8 

0.025 2 

0.0052 2 

991208 

18 2 

1.53 2 

0.056 2 

0.035 2 

991216 

4.7 2 

1.36 2 

0.014 2 

0.018 2 

000301C 

27 4 

1.43 2 

0.062 2 

0.072 2 

000418 

27 2 

2.04 2 

0.076 2 

0.0066 2 

000926 

22 4 

2 5 s 2 , 3 , 8 

0.125 2,3 

0.0435 2,3 

010222 

1.7 2 

1.695 2,8 

0.43 2 

0.000067 2 

011121 


2.5 s 


0.5 s 

011211 

3 5 

2.4 5 

0.0025 s 

0.01 s 

050801 


2.64 9 



050802 


2.62 9 



050904 



0.0309 7 


051109A 


2.08 9 



060124 


2.02 9 



060729 


2.22 9 



061121 


1.88 9 



090323 



0.07 10 

0.0089 10 

090328 


2.26 10 

0.11 10 

0.0019 10 

090902B 

0.00056 10 

2.21 10 

0.13 10 

0.33 10 

090926A 


2.13 10 

0.33 10 

0.081 10 


1 


Wiiers Galama (£l99£)I) 

2 


Panaitescu & Kumar (20021 

3 


Yost et al.| ( 20o3) 

4 


Chevalier 8z Li 1 2004|) 

5 


Piro et al. 1 20051) 

6 


van der Horst et al. < 20051) 

7 


Gou et al. ( 20071) 

8 


Starling et al.| ( 20081) 

9 


Curran et al.l d2C^^) 

10 

Cenko et al.l (201]]) 


Where multiple references are given, the value shown is the nu¬ 
merical average of the reference values. 


distribution’s median value of 0.88 is not dissimilar to the 
fiducial density of 1.0 cm -3 . 


3.5 Electron energy index p 

Previous studies have derived a larger number of p val¬ 
ues, as well as possible distributions of this parameter. Sev¬ 
eral theoretical works concluded that p should have a near- 
universal value of ~ 2 . 2 - 2.3, with very little variation 
llAchterberg et al.1 l200ll : iLemoine fe Pelletier! 120031 1 . Obser¬ 
vationally, p varies much more widely, as seen in Table |T] 

A Gaussian distribution that fits the data in Table [T] 
has /r = 2.2, a = 0.4. The K-S test of this data against a 
distribution with these parameter^] yields a 99% probability 
that this data came from such a distribution. Further, the 
mean of 2.2 is consistent with the fiducial value of 2.3. 


3 scipy.stats.kstest(data, ‘norm’, args = (2.2, 0.4)) 



Figure 4. Comparison of o bserved GRB radio afterglows (dot¬ 
ted lines) llChandra fe Frail |2012| 1 and simulated Pop III GRB 
afterglows (solid lines) at 8.5 GHz. 

3.6 Electron (e e ) and magnetic (ea) energy 
fractions 

Observationally derived values for t e are shown in Table 
□ At first glance, the values roughly conform to a slight 
variance around the fiducial value of 0.1. 

Using the K-S test, we find a 74% probability that the 
observed values of e e are from a lognormal distribution with 
e M = 0.09 and <7 = 1.2. The median value compares well 
with the fiducial e e of 0.1. 

Using the K-S test, we find a 98% probability that the 
tB values shown in Table |T] are from a lognormal distribution 
with e M = 0.012, a = 2.6. We are again assured by the 
proximity of the median to the fiducial tB of 0.01. 


4 RESULTS 

With the parameter distributions described in Section 3, 
we use the GRB afterglow model described in Section 2 to 
simulate GRBs from Pop III stellar collapse. The proba¬ 
bility distribution functions of the parameters are used to 
generate pseudo-random values, which are then fed into the 
simulation to produce lightcurves. 

Firstly, by testing those simulations where the observa¬ 
tion angle is within the initial opening angle (the GRB is 
aligned with the observer), and comparing to observed af¬ 
terglows, we determine that the afterglows of aligned Pop 
III GRBs will be observable. Using the 8.5 GHz radio af¬ 
terglow data provided by Poonam Chandra (2014, private 
communication), the comparison of simulated Pop III GRB 
afterglows to observed GRB afterglows is shown in Figure^ 
The figure shows that the simulated lightcurves are similar 
to observations, so the afterglows of Pop III GRBs should 
easily be detectable at frequencies near 8.5 GHz with exist¬ 
ing instruments. 

As shown in Figure [4] there is nothing in the shape 
of the afterglow radio lightcurves to distinguish between 
GRB progenitors. The lightcurves from our simulations pref¬ 
erentially occupy the higher peak flux density region of 
the spread of observed lightcurves, but no individual burst 
would be noticeably different. One may have expected that 
the higher energy budget presumed available to Pop III 
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Figure 5. 8.5 GHz radio afterglows of simulated Pop III GRBs 
when the observer is on-axis (black lines) and off-axis (grey lines). 

GRBs would produce brighter afterglows, and indeed in the 
local ISM frame this would be the case. However, the high 
redshift of these GRBs then reduces the observed flux den¬ 
sities to values typical of local GRBs. 

The simulations show that the radio afterglows of Pop 
III GRBs should not appear very different from normal 
GRBs. Conversely, it also means that an observed radio af¬ 
terglow could potentially be a high redshift Pop III GRB. 

Having shown that Pop III GRB afterglows are de¬ 
tectable, in theory, we also consider the contribution of un¬ 
aligned ‘orphan’ GRBs, where the observing angle is initially 
outside the jet opening angle. Figure [5] shows that the flux 
densities can reach detectable levels after the jet break . 

4.1 Detection by SKA precursors 

In this study we consider two Square Kilometre Array 
(SKA) precursor instruments; the Murchison Widefield Ar¬ 
ray (MWA) and the Australian SKA Pathfinder (ASKAP). 

MWA is a low frequency wide field phased array instru¬ 
ment, operating in the frequency range 80-300 MHz. The de¬ 
fined characteristic performance of this instrument is given 
at the central frequency of 150 MHz with a 30 MHz band¬ 
width, where it has a nominal field of view of 6J.0 de g 2 a nd 
angu lar resolution of between ~ 2 ~ 3 arcmin dTingav et al.l 
120131) . The sky noise at this frequency, combined with the 
angular resolution, results in MWA being largely confusion 
limited. The 5 ct detection limit is 6 mjy, which the instru¬ 
ment achieves in a matter of minutes. When simulating the 
MWA lightcurves we integrate over a band 30 MHz wide 
centred on 150 MHz, using 6 mjy as the sensitivity limit. 

ASKAP’s operational frequency range is 0.7-1.8 GHz, 
with a bandwidth of 300 MHz. With the instrument config¬ 
ured for an angular resolution of 18”, it has a 5a po i nt source 
sensi tivity of 29 fiJy in a 1 hour observation (IJohnston et all 
120091 ). In our simulations of the ASKAP lightcurves we inte¬ 
grate over a band 300 MHz wide centred on 1.25 GHz, and 
use 29 /rjy as the sensitivity limit. 

Figure [6] shows the peak flux densities of the simulated 
afterglows plotted against redshift for ASKAP and MWA. 
Grey points indicate that the peak flux density is less than 
the detection limit, and black points are above the detection 
limit. 


Figure 6. Peak flux densities of simulated lightcurves of Pop 
III GRB radio afterglows in the frequency bands of ASKAP and 
MWA plotted against redshift in the most optimistic case. Black 
points show afterglows exceeding the detection threshold, while 
grey points are undetectable afterglows. 

Figure[6]shows that ASKAP will detect many more Pop 
III GRB afterglows than MWA. 34.8% of all afterglows could 
be detectable by ASKAP. Only 4 of the three thousand sim¬ 
ulated afterglows achieved a peak flux density greater than 
6mJy at 150MHz, a detectability rate of 0.13%. Closer anal¬ 
ysis of the data reveals 89.4% of aligned afterglows will ex¬ 
ceed ASKAP’s detection limit, and that the majority (70%) 
of those will do so within one day after the burst. 

The simulated lightcurves whose peak flux are shown in 
Figure [6] assume a very optimistic scenario for Pop III SFR 
and IMF. Properly investigating the possibility of Pop III 
GRB radio afterglow detection requires that we also consider 
the worst case (and likely more realistic) sce nario. The most 
pessim istic derived Pop III SFR is that of iKulkarni et al.1 
(l2013h . where the redshift of peak Pop III star formation 
occurs at higher z, mostly at z > 10, and never at 2 < 7.5. 

The Kulkarni SFR also assumes the least top-heavy, 
thus for us the most pessimistic, Pop III IMF of a Salpeter 
power law distribution with a low mass cut-off of IMq. The 
pessimistic SFR and IMF yield new probability distribution 
functions for redshift and Ei so . The simulated afterglows’ 
peak flux densities with these PDFs are shown in Figure 0 

The figure shows that even in the pessimistic scenario 
of Pop III formation, the resultant GRB radio afterglows 
will be detectable by ASKAP. There is a reduction in the 
proportion of detectable afterglows; in this scenario, 82.7% 
of aligned afterglows and 26.9% of all afterglows exceed the 
limit of ASKAP. Only a single simulation (0.03% of all af¬ 
terglows) exceeded the detection limit for MWA. 


4.2 Pop III GRB afterglow detection rates 

Our simulations give a representation of the detectability of 
Pop III GRBs. To assess the probability of making a de¬ 
tection, we must estimate the intrinsic rate at which these 
GRBs occurs. I 11 constructing our sets of simulated GRB 
afterglows, we have used theoretical SFRs and IMFs, both 
of which affect the formation rate of GRBs. Using the ap¬ 
proach of Bromm & Loeb (2002), the GRB rate above a 
certain redshift is as follows: 
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Figure 7. Peak flux densities of simulated lightcurves of Pop III 
GRB radio afterglows in the bands observable by ASKAP and 
MWA plotted against redshift in the most pessimistic scenario 
of Pop III SFR and IMF. Black points show afterglows exceed¬ 
ing the detection threshold, while grey points show undetectable 
afterglows. 



Ngrb(> z) 


47TcA£ obs r°° _ r/GRB 

Ho J z (1 + z) 3 + z ) 3 + Qa 


(32) 

where >]/» is the Pop III SFR, Af 0 b s is the observation time 
interval, set to 1 yr to cancel the yr~ 2 component of the SFR, 
Ho is the Hubble constant, and r/ grb is the GRB formation 
efficiency: 


IJGRB = /GRB 


Sm U l <Krn)r(m)dm 
Sm” mcj>{m)dm 


(33) 


where My and My are the upper and lower lim its of <j>(m ), 
the IM F. r(m) is the probability of satisfying the lYoon et al.l 
(l2012h rotation speed criteria as a function of mass, /grb 
is the GRB formation fraction, being the fraction of poten- 
tial progenitors which ac tually produce GRBs, set at 0.001 
dLanger fe Normanll2006l 'l. 

In the most optimistic case of Pop III GRB formation, 
we us e the most optimistic Pop III SFR of lde Souza et ahl 
J201lh . combined w ith a Gaussian IMF wit h /x = 55M© 
and a = 15Mm (IScannapieco et al.l 120031 ') as used in 
Ide Souza et al.l d201lh . With this IMF. r/GRB is 6.35 x 
10 _6 Mq . Applying this and the optimistic SFR to equa¬ 
tion [321 for all redshifts results in an intrinsic Pop III GRB 
rate of 4.55 x 10 3 yr _1 . We can scale this rate by the fraction 
of the sky covered by the ASKAP, and find that Pop III 
GRBs would occur at a rate of 3.31 yr -1 within its 30 sq. 
deg. FOV. 

Multiplying the average rate of GRBs per day per FOV 
by the average number of days that the afterglows are de¬ 
tectable, results in the average number of detections per 
image. Where this number is greater than 1, the implication 
is that, at any given time, there should be multiple simulta¬ 
neously observable GRB afterglows. Where this number is 
less than 1, it implies multiple observations will be required 
per detection. 

We only perform this analysis for ASKAP, for as shown 


KT 1 



Observer Time (days) 


KT 1 



Observer Time (days) 


Figure 8. Simulated lightcurves of Pop III GRB radio afterglows 
at a frequency band 1.10-1.40 GHz, observable by ASKAP. Grey 
segments of the lightcurves are below ASKAP detection thresh¬ 
old of 29/i.Jy; black segments are above this threshold. Top The 
optimistic scenario. Bottom The pessimistic scenario. 


in Figures [S] and [3 the probability of MWA detecting a Pop 
III GRB afterglow is negligible. 

Figure [8] (Top) shows the simulated GRB afterglows re¬ 
sulting from using this SFR and IMF. Analysing the data, 
we find that the average length of time that a Pop III GRB 
afterglow remains detectable by ASKAP is 6.60xl0 5 days, 
inclusive of those which never exceed the detection thresh¬ 
old. If we exclude the invisible GRBs, then the remainder 
are detectable by ASKAP for an average of 1.90x10® days. 

In this optimistic case, for ASKAP we have 9.06xl0~ 3 
Pop III GRBs per day multiplied by 6.60 x 10 5 days, resulting 
in a detection number of 5.98xlO 3 . These ~ 6000 Pop III 
GRB afterglows in every 30 sq. deg. image taken by ASKAP 
would, on average, be separated by 4.8 arcmin. This scenario 
should therefore be easily proven or disproven by ASKAP. 

The most pessimistic case of Pop III GRB f ormation 
rate c onsists of the Pop III SFR derived bv iKulkarni et ahl 
(120131) , adopting a Salpeter power law IMF in the range 
1-100 Mq. This IMF combined with the aforementioned for¬ 
mation criteria results in t/grb = 1-82 x 10 _6 Mg. Applying 
this efficiency and the pessimistic SFR to eauation l32l results 
in an intrinsic Pop III GRB rate for all redshifts of 7.95 yM 1 . 
This equates to a rate of 1.58x 10 _5 yr _1 within ASKAP’s 
FOV. 

From Figure [8] (Bottom), the average time an afterglow 
can be detected by ASKAP is 6.72xl0 5 days. The average 
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Figure 9. Observed GRBs E; so as a function of redshift, roughly 
d 2 

following the expected yp- relation. 


number of detectable Pop III GRB afterglows per ASKAP 
image is 10.6, with an average angular separation between 
afterglows of 1.2 degrees. 

Since the proportions of afterglows detectable by MWA 
is negligible in both pessimistic and optimistic scenarios, we 
do not calculate detection rates for this instrument. 


4.3 Detection of prompt gamma-ray emission 


So far we have not been concerned with the orientation of 
the GRBs, treating the small subset of directed bursts right 
along with those undirected bursts which only become vis¬ 
ible at late times. In terms of performing a blind search, 
whether or not the observer was within the initial GRB 
opening angle is irrelevant. We now look at the prospects 
of a triggered follow up of Pop III GRBs from the prompt 
gamma-ray detection by Swift. 

The first question is whether or not the prompt gamma- 
ray emission from a Pop III GRB will be sufficiently lu¬ 
minous to trigger a GRB satellite. We assume that peak 
gamma-ray luminosity scales with Ei SO and we observe that 
the minimum detectable peak luminosity, and therefore min¬ 


imum detectable Ei so , increases with redshift as 


di 


1 + z 


• Fig¬ 


ure [9] plots the limiting Ei so with redshift following the ex- 
d 2 

pected Ei so lim oc ——. To ensure a consistent sample of 
1 + z 

Eiso data, we use the Butler online catalogue Swift BAT In- 
te grated Spectral Par a meter^- This catalo gue, an extension 
of lButler at alJ (l2007lf : feutler et al.l (l201Cll l. circumvents the 
nominal BAT upper energy of 150 keV to produce values of 
Ei ao through a Bayesian approach. 

The figure shows that at z = 35, which is the upper limit 
for the theoretical Pop III SFRs adopted in this paper, the 
minimum observable Ei so is less than 10 54 ergs. It is less than 
the minimum Ei SO of Pop III GRBs assumed in this paper 
(see Figure [2]), so in this scenario the prompt emissions will 
be sufficiently luminous to trigger Swift. 

The second question is the rate at which both Swift is 
within the initial opening angle of a Pop III GRB, and that 


4 http://butler.lab.asu.edu/Swift/bat_spec_table.html 


GRB is within the FOV of the satellite. Within Swift's 2 str 
FOV, Pop III GRBs occur at a rate of 724 yr _1 in the opti¬ 
mistic case, and 1.27 yr _1 in the pessimistic case. Through 
our simulations, the observer was within the initial opening 
angle of approximately 4% of all GRBs. In the optimistic 
case Swift would detect 29 Pop III GRBs yr _1 , and in the 
pessimistic case 0.05 yr _1 (one Pop III GRB every 20 years). 

Following Swift detection, ASKAP would be capable 
of detecting ~90% of the 29 GRBs yr _1 (26 yr _1 ), in the 
optimistic case, and have an ~83% chance of detecting a 
1/20 year Pop III GRB in the pessimistic case. MWA should 
not be considered to follow up Pop III GRB prompts in 
any case, as our simulations have shown that a) the chances 
of achieving sufficient flux density for MWA to detect are 
extremely low, and b) the sufficient flux density occurs a 
very long time after the burst. 


4.4 Optimal observing cadence for detection 

The optimal cadence to search for slow transients such as 
GRB radio afterglows are those that have the highest aver¬ 
age flux variance between images. We consider an observa¬ 
tion plan lasting ten years, and populate a number of theo¬ 
retical FOVs with simulated afterglow lightcurves. We took 
a range of observation intervals, from observing every day 
for the full ten years, to having only two observations sep¬ 
arated by ten years. Over this range we find the average 
changes in afterglow flux density over a single time interval. 

For the optimistic case with ASKAP, we find that the 
average change in flux density increases with the length of 
the interval, up to the maximum interval of 10 years, where 
the average change in flux density was ~ 6 x 10~ 5 Jy. The 
pattern continued in the pessimistic case, where the maxi¬ 
mum average flux change of ~ 2.7 x 1CU 5 Jy occurs at the 
maximum interval of 10 years. 

We find the optimum cadence for these searches with 
instruments like ASKAP will be greater than ten years. For 
practical purposes, the results indicate that when it comes to 
a blind search for Pop III GRB radio afterglows, one should 
look at the longest feasible cadence. 


5 CONCLUSIONS 

We find the radio afterglows of Pop III GRBs are detectable 
by ASKAP, via the application of a realistic afterglow model, 
under the following assumptions: 1: Pop III stars formed at 
high redshift, with higher masses and therefore shorter lives 
than Pop I/II stars. 2: Pop III stars could produce GRBs in 
the manner of stellar core collapse, providing certain criteria 
of mass and rotation were met. 3: The resulting GRBs would 
be more energetic than those from Pop I/II stellar collapse. 
4: The microphysical parameters affecting GRBs and their 
afterglows do not change with stellar population. 

Our simulations show that the radio afterglow 
lightcurves of high redshift Pop III GRBs will be indistin¬ 
guishable from those of low redshift Pop I/II GRBs. This 
undermines the primary motivation for searching for Pop 
III GRBs: the study of Pop III star formation history. The 
only indication that a radio afterglow might belong to a Pop 
III star would be if the GRB were at high redshift (z > 5). 
The fact that the radio emission is similar to that of a low 
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redshift GRB (z < 5) would indicate Pop Ill-type energy. 
The distinction cannot be based solely on the radio after¬ 
glow lightcurves, as unfortunately there is no way to use 
radio afterglow data to determine redshift. 

Even if radio observations could identify high redshift 
and orphan GRB afterglows, those observations would not 
by themselves provide evidence that the progenitor were 
Pop I/II or III. Observations at higher frequencies will be 
necessary, but at these redshifts the Gunn-Peterson trough 
of Lyman-a absorption rules out frequencies higher than 
near-infrared. With advanced space infrared telescopes (e.g. 
James Webb Space Telescope) it may be possible to obtain 
spectra of high redshift GRB afterglows, which would reveal 
the composition of their environments. However, given that 
optical/infrared afterglows tend to decay quickly, it is likely 
that such observations would only be possible for GRBs di¬ 
rected at and triggering a GRB satellite. 

It is expected that the peak gamma-ray luminosity of 
Pop III GRBs will be high enough to trigger a GRB satel¬ 
lite, such as Swift, despite their higher redshifts. The rate 
of prompt detection relies on the alignment rate scaled by 
the sky fraction covered by GRB satellites. Swift could de¬ 
tect as many as 29 Pop III GRB prompts yr _1 in the most 
optimistic case, and as few as one every 20 years in the pes¬ 
simistic case. 

Simulations of the afterglow lightcurves of Pop III 
GRBs in the optimistic case show that 89% of the aligned 
bursts will be detectable by ASKAP, 70% of which within 
1 day. ASKAP’s sensitivity lends it to timely follow-up of 
GRB prompts, although rapid follow-up interrupting its nor¬ 
mal observation schedule should not be necessary due to the 
longevity of the afterglows. Observations made some tens of 
days after the GRB may identify an afterglow with a Pop 
III progenitor if it has a brightness typical of a low redshift 
Pop I/II progenitor and is shown to be high redshift (z > 
5). 

Our simulations do not find any afterglows of aligned 
Pop III GRBs detectable to MWA. The total number of 
detections was negligible (0.13%), and only possible when 
near their peak flux density, which is extremely late at MWA 
frequencies. MWA is unsuited for follow up of Pop III GRB 
afterglows, or indeed GRBs in general. 

Altogether 35% of all aligned and unaligned Pop III 
GRB afterglows will be detectable by ASKAP in the op¬ 
timistic case. Since, on average, only 4% of all GRBs will 
be aligned with the observer, the feasibility of detecting the 
unaligned ‘orphan’ Pop III GRB afterglows with these in¬ 
struments is important. The average number of detections 
per image is the product of the intrinsic rate of Pop III 
GRBs within the FOV with the average time an afterglow 
is detectable. For ASKAP this is 5.98xl0 3 detectable af¬ 
terglows per image, a number which is unrealistically high. 
The number of detectable afterglows predicted in this case is 
high enough to be relatively easily disproved, and we expect 
ASKAP transient searches to dismiss this upper limit. 

In the pessimistic case, analysis of the simulations shows 
ASKAP capable of detecting 83% of aligned afterglows and 
27% of all afterglows, with an average of 10.6 detections per 
image. MWA detections are negligible (0.03%). 

The average change in afterglow flux density over a time 
interval increases with the length of that interval, at least up 
to the maximum tested interval of ten years. This average 


flux variance seems to be of the order of the flux density 
limit of detection. In performing a blind or archival search 
for orphan afterglows, one should use the longest feasible 
cadences and search for flux decay. 

While detecting orphan radio afterglows may be possi¬ 
ble, our results show that it will not be possible to distin¬ 
guish progenitors from orphan radio afterglow data alone. 
Combining such data with other detections may, in principle, 
reveal new insight into progenitors. For instance, detections 
of an orphan afterglow in multiple frequencies from sub-mm 
to radio, and a non-detection in optical/infrared frequencies 
may be due to redshifted Lyman-a absorption indicative of 
a high redshift and high energy GRB progenitor. The key 
question would be the proportion of Pop III GRBs whose or¬ 
phan optical/infrared afterglows would be detectable if not 
for redshifted Lyman-a absorption. This question we will 
investigate in a follow-up study. 
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